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A non-zero CP-violating 6 parameter is treated in the domain model which assumes a cluster-like 
vacuum structure whose units are characterised in particular by a topological charge which is not 
necessarily an integer number. In the present paper we restrict consideration to rational values 
of the charge. The model has previously been shown to manifest confinement, spontaneous chiral 
symmetry breaking and the absence of an axial U{1) Goldstone boson. We find that the specific 
structure of the minima of the free energy density of the domain ensemble forces a 27r-periodicity 
of observables in 9 for any number of light quarks, that vacuum doubling occurs ai 6 = n for any 
Nf > 1 and any value of topological charge q. These features are in agreement with expectations 
based on anomalous Ward identities and large Nc effective theories. We find also additional values 
of depending on q for which vacuum doubling occurs. 
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I. INTRODUCTION 

Explicit CP violation can be introduced in quantum 
chromodynamics by the inclusion of the so-called theta 
term in the action. In Euclidean space this amounts to 

2 p 

Se = ^q9, 9 = ^ y d'^xF^^F^'^. (1) 

A remarkable feature of such a term in the presence of 
spontaneously broken chiral symmetry is the Dashen phe- 
nomenon iJj: this explicit breaking of CP can become 
spontaneous due to vacuum doubling at the point 9 = n. 
This value falls outside the physically relevant range 
< 9 < 10~^ to which 9 is constrained by the neutron 
dipole moment and current algebra 0. Nonetheless, 
the Dashen phenomenon is a constraint on self-consistent 
models of the non-perturbative QCD vacuum. The state- 
ment of the problem related to the Dashen phenomenon 
as well as the question about periodic dependence of ob- 
servables in 9 acquire full sense only if the values of q are 
not restricted to integers. In this paper we use the do- 
main model jJI of the QCD vacuum to show how specific 
properties of quark field configurations summed up in the 
partition function can lead to the periodic dependence on 
9 and to the Dashen phenomenon for any rational values 
of q. 

The approaches which originally demonstrated the 
Dashen phenomenon in the context of QCD are those 
of anomalous Ward identities and effective chiral La- 
grangians 1^ in the limit of large number of colours 
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Nc, which we shall describe below in more detail. Sub- 
sequent discussions of the Dashen phenomenon include 
II 11 IS ES El- This phenomenon for 6i is inti- 
mately related to both the mechanism of spontaneous 
chiral symmetry breaking and the non-appearance of an 
axial U{1) pseudoscalar Goldstone boson as shown by 
the original works revealing it in QCD. In the key 
tool is the generalisation to the axial U{1) channel of 
Ward identities in which the true vacuum of the the- 
ory is unknown, but for which divergences of current- 
current expectation values in this vacuum can be related 
to hadron spectroscopy via vacuum symmetry proper- 
ties. For example, the light pion mass as input into 
the fiavour SU{2)l x SU{2)ji Ward identity (which gives 
the Gell-Mann-Oakes-Renner relation) leads to a non- 
zero chiral condensate as the output. In the axial U{\) 
channel (where the anomaly figures) and = the phe- 
nomenological input is now the absence of a light meson, 
with the output that topological charge q must be frac- 
tional. With 6* 7^ the corresponding output starting 
with the same phenomenological input is 27r-periodicity 
in 9 and the property of vacuum doubling and sponta- 
neous CP breaking at = tt. In the large Nc approach 
the same data is turned around. An effective Lagrangian 
for mesons is written down for Nc oo which assumes 
chiral symmetry breaking and the anomaly. The out- 
put for 9 = Q includes the famous relation between the 
topological susceptibility and various meson masses. For 
9 ^ Q again the Dashen phenomenon emerges. 

However these approaches do not unveil the actual 
mechanism of non-perturbative vacuum properties such 
as confinement and spontaneous chiral symmetry break- 
ing. The domain model |J| is a bottom up approach to 
these features: a particular vacuum structure is intro- 
duced into the formalism explicitly, here based on a sta- 
tistical ensemble of domain-like gluon fields, and out of 
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this both vacuum and mesonic properties are derived. In 
previous works both confinement Q and chiral symme- 
try reahsation 0, 0| have been studied for this vacuum 
ansatz. The anomaly contribution to the free energy sup- 
presses continuous axial U{1) degeneracy in the ground 
state, leaving only a discrete residual axial symmetry. 
This discrete symmetry and flavour SU{Nf)L x SU{Nf)R 
chiral symmetry in turn are spontaneously broken with 
a quark condensate arising due to the asymmetry of the 
spectrum of Dirac operator in the domain background. 
An estimate of pseudoscalar and vector meson masses 
has been computed showing the typical pattern of spon- 
taneous breakdown of SU{Nf)L x SU{Nf)ji chiral sym- 
metry with the correct splitting between tt, rj and rj' 
mesons In this paper we test the model against 

the above predictions of spontaneous CP breaking and 
generally explore its features for non-zero 9 parameter. 
The Dashcn phenomenon is studied by directly comput- 
ing the 0-depcndendence of the free energy density and 
chiral condensate within the domain model. 

The results of this paper are that indeed the vacuum 
doubling at 9 = tt occurs in the domain model for any ra- 
tional values of the topological charge q and any number 
of flavours. Thus the domain model provides an explicit 
example of simultaneous realisation of confinement and 
chiral symmetry breaking where the Dashcn phenomenon 
is manifest. Within the model, the topological charge is 
not restricted to any ad hoc value, and in the general case 
there are critical values in addition to 6* = tt, which disap- 
pear when q is set to 1 /Nf (which is the value of q emerg- 
ing in jSj] for the case of Nf degenerate quark flavours). In 
particular, in 0| the domain model parameters were fixed 
from the string tension giving q = 0.15. For this topolog- 
ical charge, the critical values 9 = — 7r/3, 7r/3, 7r(mod27r) 
appear for which the vacuum doubling occur. As with 
vacuum degeneracy for discrete 9 and 27r-periodicity 
of observables such as the condensate arises as output 
rather than being assumed. Unlike 0,0, these results 
emerge without the absence of a massless U{1) boson 
being input. 

We first briefly review the domain model for 9 — and 
thereafter examine the realisation of chiral symmetry for 
9 0. In the appendices we review and compare our 
results with the approaches of and . 



partition function for iV — > oo domains of radius R 



II. REVIEW OF CHIRAL SYMMETRY IN THE 
DOMAIN MODEL 



Z = 



/OO ^ /> /> 



xe 



(2) 



where the functional spaces of integration J-q and are 
specified by the boundary conditions {x — Zif' — R^ 



[x) = 0, 



i^,{x)e"'^''^^\x) = i^^'^x), 



(3) 



(i) iQ75 



ifli{x) 



ij^'^ix). 



Here fii = n^t" with the generators f'' of 5'f7c(3) in the 
adjoint representation and a is a random chiral angle 
associated with the chiral symmetry violating boundary 
condition Eq.(|3Jl in the presence of one species of quark. 
The generalisation to Nf flavours is given in the next 
section. The thermodynamic limit assumes V, N ^ oo 
but with the density = N/V taken fixed and finite. 
The partition function is formulated in a background field 
gauge with respect to the domain mean field, which is ap- 
proximated inside and on the boundaries of the domains 
by a covariantly constant (anti-)self-dual gluon field with 
the field-strength tensor of the form 



N 



ix~z,f/R^), 



with bI/JbI'^ 



B'^Sjyp- Here z^ are the positions of the 



centres of domains in Euclidean space. 

In order to reflect continuity of colorless quark cur- 
rents at the positions of pure gauge singularities in the 
picture of quark and gluon configurations, which is ap- 
proximated by the domain model, the chiral angle a is 
taken to be the same for all domains in the represen- 
tation Eq. This is a refinement of the formulation 
presented in (4illlll^. The measure of integration over 
parameters characterising domains is 



27r 

1 f d^Zi f 
A8tt Jy V J Jo 
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(4) 



For motivation and a detailed description of the do- 
main model we refer the reader to 0. The essential def- 
inition of the model is given in terms of the following 



where {9i, ipi) are the spherical angles of the chromomag- 
netic field, uji is the angle between chromoelectric and 
chromomagnetic fields and is an angle parametrising 
the colour orientation. 
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This partition function describes a statistical system of 
domain-Hke structures of density , where the volume 
of a domain is w = Each domain is characterised 

by a set of internal parameters and with internal dy- 
namics represented by fluctuation gluon Q^*-* and quark 
-^f^'^ fields. It respects the symmetries of the QCD La- 
grangian, since the statistical ensemble is invariant under 
space-time, colour gauge and particularly chiral symme- 
tries. Thus the model involves two free parameters: the 
mean field strength B and the mean domain radius R. 
These dimensionful parameters break the scale invariance 
present originally in the QCD Lagrangian. In principle, 
they should be related to the trace anomaly of the energy- 
momentum tensor il9, 20] and, eventually, to the funda- 
mental scale Aqcd- Knowledge of the full quantum ef- 
fective action of QCD would be required for establishing 
a relation of this kind. Within this framework the gluon 
condensate to lowest order in fluctuations is AB^ and the 
topological charge per domain is g = B^R^ /16. 

An area law is obtained for static quarks. The reason 
for this is the finite range of gluon correlations implicit 
in the model which will figure in all the phenomena we 
consider. Computation of the Wilson loop for a circu- 
lar contour of a large radius L ^ R gives a string ten- 
sion a = BfiirBR^) where / is given for colour SU{2) 
and SU{3) in ||. Estimations of the values of these 
quantities are known from lattice calculation or phe- 
nomenological approaches and can be used to fit B and 
R. As described in ^ these parameters are fixed to be 
Vb = 947MeV,i? = (760MeV)-i = 0.26 fm with the 
average absolute value of topological charge per domain 
turning out to be q = 0.15 and the density of domains 
v~'^ = 42fm~''. The topological susceptibility then turns 
out to be X = (197MeV)^, comparable to the Wittcn- 
Veneziano value 6]. The eigenvalue problem 
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with boundary condition Eq. © was studied in l3\ . With 
the domain background field the solution to this problem 
reveals an asymmetric spectrum, exhibiting the broken 
chiral symmetry through the bag-like boundary condi- 
tion, and thus none of the eigenmodes is chiral. However 
at the centre of domains all modes are chiral and the 
sign of their chirality depends on whether the underly- 
ing gauge field is self-dual or anti-self-dual. In |l2| we 
computed the distribution of values of the local chiral- 
ity parameter of [l7| in a chirally symmetric ensemble 
revealing qualitatively similar behaviour to the double- 
peak structure seen on the lattice Ts], which is taken to 
be indicative of spontaneously broken chiral symmetry. 

The a-dependent part of the free energy density was 
computed in 'l3l| using zeta function regularisation, with 
an imaginary part arising 



— ±— arctan(tan(a)) 

V 



(6) 



where q is the absolute value of topological charge in a 
domain, and overall sign (— corresponds to an (anti- 



)self-dual domain. The charge is not integer here in gen- 
eral but the anomalous term is nonetheless im periodic 
in a. This is the Abelian anomaly as observed within 
the context of bag-like boundary conditions by and 
is consistent with [l5| albeit not generated from purely 
chiral zero modes but the chiral properties of non-zero 
modes. 

The part of the free energy density J-' relevant for the 
present consideration of an ensemble of A'^ ^ oo domains 
with both self-dual and anti-self-dual configurations takes 
the form 

/oo 
da [cosvQF{a)f 
-OO 

= ^f exp(A^ln[cos(^;5F(a))]) + 0(1/Af). 

Here summation goes over the infinite set amin of de- 
generate minima of the free energy density, which are 
achieved at 

v^F{ajnin) — O(mod 2n). 

In the thermodynamic limit N ^ oo the solutions to 
the above equation for amin correspond to a degenerate 
set of vacua connected by discrete chiral transformations, 
which will be given in the next section. 

Thus for massless quarks, a discrete subgroup of [/^^(l), 
rather than the continuous C/a(1) itself, represents a sym- 
metry of the vacua. The anomaly defines those chiral an- 
gles which minimise the free energy so that the full J7a(1) 
group is no longer reflected in the vacuum degeneracy. It 
should be stressed here that this residual discrete degen- 
eracy ensures a zero value for the quark condensate in 
the absence of mass term or some other external chiral- 
ity violating sources. 

In the presence of an infinitesimally small quark mass 
the a dependent part of the free energy of a self-dual 
domain is modified by the term linear in mass and takes 
the form (for details see 13]) 



F 



i-arctanftanfa)] — mHe'" + 0(ni ), 

V 



E 



7r2i?3 ^ fc + 1 

k—l,Z — Zi ,Zi,Z2 



[M{l,k + 2,z) 



k + 2 



M{l,k + 3,-z) - 1 



(7) 



where the quantity H > comes from the spectral asym- 
metry term 77(1), 



77(s)=^sgn(A)/]A]^ 



appearing due to the asymmetry of the spectrum [l6l] . 
The free energy of an anti-self-dual domain is obtained 
via complex conjugation. 

This discrete chiral symmetry of the massless case is 
spontaneously broken and switching on the quark mass 
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FIG. 1: The scalar quark condensate as a function of the 
parameter 6 for Nf = 1 and g = .15 in units of K. The dashed 
hues correspond to discrete minima of the free energy density 
which are degenerate for m = 0. The sohd hne denotes the 
minimum which becomes global for a given 6 in the presence 
of infinitesimally small mass term. 



selects one of the vacua. For our conventions of boundary 
condition and mass term the selected minimum is at a = 
0. It is important that values for the angles Ofmin are 
not modified by the leading order term, linear in m. The 
quark condensate is extracted from the free energy in the 
standard way 



— lim lim — — ■— — e 

rn^O N^oo vN dm 



vNT{m) 



= -H 



(8) 



and takes the value {%p{x)%lj{x)) = -(237.8 MeV)^ for 
the values of field strength B and domain radius R fixed 
earlier by consideration of the pure gluonic characteristics 
of the vacuum - string tension, topological succeptibility 
and gluon condensate. 

For Nf > 1 quark flavours the fermion boundary con- 
dition in Eq.Q, explicitly breaks all chiral symmetries, 
flavour singlet and non-singlet (see also Thus in- 

tegrating over all a does not suffice to provide for the 
full chiral invariance of the ensemble of quark configu- 
rations contributing to the partition function. Rather, 
the boundary condition must be generalised to include 
flavour non-singlet angles, a a -I- /3T°/2, with T" 
the N'j — 1 generators of SU{Nf). Then integration over 
N'j angles a and /3° must be performed for a fully chiral 
symmetric ensemble. The spectrum of the Dirac opera- 
tor can be found now quite analogously to the one-flavour 
case, except that the boundary condition mixes flavour 
components and an additional projection into flavour sec- 
tors is required in order to solve Eq.lQ. The calculation 
from this point will be repeated for non-zero 9 in the 
next section, but it suffices to summarise here the re- 
sult for 9 = emerging in jlSi]: the Abelian nature of 
the anomaly meant that the non-singlet chiral angles im- 
plicit in this procedure drop out of the ensemble free en- 
ergy. Thus there are only Nj — 1 continuous directions in 



the vacuum. The expectation that there should be onl; 
N'j — 1 pseudo-Goldstone bosons has been verified in 
for Nf = 3 by an estimation of the meson spectrum. 
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FIG. 2: The scalar quark condensate as a function of Q for 
N j = 3 and g = .15 in units of K. The meaning of dashed and 
solid lines is the same as in Fig. Qualitatively the picture 
is the same as for N j = 1: the same periodicity and critical 
values of 5, which is achieved by increased number of different 
discrete minima. 



III. FREE ENERGY IN THE PRESENCE OF 
NON-ZERO Q 

Now we include the CP-violating parameter in the 
model by the additional term Eq. JQ) in the action, which 
contributes a pure phase to the weight factor in Euclidean 
space. Integrating over N f fermions with infinitesimally 
small masses mi = • . • = jntq = m in a domain ensemble 
gives for the free energy density per domain 



F = —V Mn cos q[WNf ~ 9] — mH cos $i, 
where 



4=1 



W. 



Bi 
Bi 



Nf 

arctan(tan$i) 

i=l 

a + B, 

for Nf = 1, 

1/31 „ _ 



for Nf = 2 



and 



Bi 



b^ + b^/Vs, B2 = -b^ + b^/V3, 
-2b^/V3, for Nf = 3. 



(9) 

(10) 

(11) 
(12) 

(13) 
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Here a is the C/(l) chiral angle, /3 are the non-singlet chi- 
ral angles for Nf—2, and 63 and are certain functions 
of the eight non-singlet chiral angles /?" for Nf^3. For 
Nf = 1, Ea. (|ll|l is just the result discussed in the pre- 
vious section. For Nf = 2 the quantities Bi and B2 in 
Ea. (|12|l arise from the projection of the quark boundary 
condition into SU{2) flavour sectors. For Nf > 2 the 
analogous functions are i?i , . . . , B^^ . For any number of 
flavours Nf the functions Bi have the property that 

Nf 

^i?, = (14) 

i=l 

which manifests the tracelessness of the flavour genera- 
tors in any basis. Thus 

Nf 

J2'^, = Nfa. (15) 
1=1 

The "arctan tan" structure in Eq. H1U|) manifests the peri- 
odicity of the free energy for arbitrary topological charge 
q. We could write 

N 

WNf = ^ arctan[tan <i>j] = Nfa (mod tt) (16) 

i=l 

which is independent of the /?° . This is the Abelian prop- 
erty of the anomaly leading to the expectation of iV| — 1 
Goldstone bosons. 
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FIG. 3: The pseudo-scalar quark condensate as a function of 
9 for Nf = 1 and g = .15 in units of N. The meaning of 
dashed and solid lines is the same as in Fig. Q 

The central effect is the chiral abelian anomaly contri- 
bution denoted as qWNf under the cosine in Eq.©. It 
should be noted that 

qWNf = Nfqa + 0{a^) for a ^ 0, (17) 

which coincides with the standard form of axial anomaly 
in the Fujikawa derivation (Note that a here is twice 
the angle of chiral transformation used in the Fujikawa's 
calculation of anomaly p^). 



The mass term is introduced as a small perturbation 
violating explicitly the chiral symmetry of the system. 
In the absence of the mass term, the minima of the free 
energy are determined by the solutions 

2t:1 Trfc „ , , , 

^kiiO)^ — + — + — (keZJeZ), 18 

Nf qNf Nf 
to the equation 

cos{qWNf - qO) = 1. (19) 

Evidently there are multiple solutions arising from the 
various periodic functions appearing in Eq.®. The in- 
dex k in the solutions reflects periodicity of tan while I 
corresponds to the periodicity of cos in Eq. ((T^ . Note 
that these solutions do not depend on the flavour non- 
singlet chiral angles, as discussed above, and thus the 
free energy displays continuous degeneracy with repect 
to these angles. 
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FIG. 4: The pseudo-scalar quark condensate as a function of 
6 for Nf — 3 and g = .15 in units of N. The meaning of 
dashed and solid lines is the same as in Fig. |2| 

The solutions Eq. H18|) give an infinite discrete set of 
degenerate minima of the free energy density. For given 
Nf and q a finite subset of these minima can be extracted 
such that all other minima are 27r-periodic copies of one 
of the already identified vacua. For instance, for Nf = 1 
and q = .15 the six different vacua correspond to all 
combinations of fc = 0, 1 and / = 0, 1, 2, while for Nf = 
2 set of "different" minima is given by combinations of 
fc = 0, 1, 2, 3 and / = 0, 1, 2. It is clear that the system of 
minima is 27r-periodic with respect to 9: 

aki{9 + 27rm)^ak'iie), fc' = fc + 2m, (20) 

and the shift in 6 can be undone by re-enumerating the in- 
finite set of minima. This periodicity does not depend on 
the value of q at all since it is ensured by re-enumerating 
of index fc without use of I. All minima are degenerate 
(the free energy is equal to zero) for zero quark mass. 



6 



Switching on the infinitesimal mass term in Eq. @ hfts 
the degeneracy due to different values of the symmetry 
breaking term 

K.. = -mH^cos($f(0)) (21) 

i 

in the free energy for different <i>f^' = aki{0) + Bi{P°-). 
This introduces non-singlct angle dependence into the 
free energy, which must now also be minimised with re- 
spect to the To this end we expand Vm in small 
fluctuations in the Bi and look for minima, bearing in 
mind the constraint Ea. (|14|l . The condition that Vm have 
stationary points leads to the condition 

sin <i>i independent of i (22) 

while the condition that this leads to a minimum forces 

cos$, >0. (23) 

Ea. (|22|l leads to ~ a+Bi being independent of i which 
means Bi = (mod 2n). On the other hand, Ea. lt^ 
is fulfilled by restricting — 7r/2 < $i (mod 27r) < tt/2. 
However, we observe that for these values of 3>i 

arctan(tan <i>i) = (mod 27r). 

Thus to guarantee a minimum with respect to non-singlet 
angles, we must have 

Nf 

Wnj- — ^ arctan(tan$i) 

i=l 

= Nf{a mod 2tt). 

Comparing to Ea. H16|l . we see that of the vacua of the 
massless free energy, aki{0), only those with k an even 
integer correspond to true minima with the mass term 
switched on. Note that these arguments only apply for 
Nf > 1. For the single flavour case the available values 
are fc = 0, 1. A similar consideration shows that for ra- 
tional q — ni/n2 the minima selected by the mass term 
correspond to I any integer if ni is odd but I even if ni 
is even. 

The result is that the 9 dependence of the quark con- 
densate is given by 

{tp{x)'ip{x))c,^,(9) = -Nf^cos{aki{0)) 

for k even and / any integer except ii q ~ ni/n2 with 
ni even for which then I is also even. This is plotted in 
Figs. Hand 121 for Nf = 1 and Nf = 3 respectively, where 
the cosine in Eq. (|21() is plotted as a function of 9 for the 
minima which are not 27r-equivalent. Critical values of 
9 correspond to the points where two minima become 
degenerate. 

The two degenerate vacua at such critical points are 
distinguished by their CP properties, which can be seen 



from the behaviour of pseudoscalar condensate as a func- 
tion of 9 

(?/;(x)75V'(x))„^,(e) = -Nfi\sin{aki{9)), 

as derived for the domain model in Appendix A. As ex- 
pected, the scalar (see Eq. (|21|l ') and pseudo-scalar con- 
densates depend on 9 through cosine and sine of aki{9) 
respectively. 

The plots for one and three flavours respectively are 
given in Figs. |3| and 01 again with q = 0.15. The pseu- 
doscalar condensate is discontinuous at the critical values 
of 9 and takes values opposite in sign for the two degen- 
erate minima: parity is thus spontaneously broken. 

In other words, we can see that for most values of theta 
the mass term selects aunique minimum of the free en- 
ergy. However, there are critical values of 9 where two 
different minima are degenerate, thus displaying a two- 
fold degeneracy of the vacuum in the presence of a mass. 

There are two conditions for critical 9. The first one is 
obviously 

aoo(^'crit) = -aki{9cTit) (mod 27r) 

where k and I should not be equal to zero simultaneously 
and without loss of generality we have taken aoo on the 
RHS. Thus 

Itt 

6»^^jj ^TYk + — (mod 27r). 

For a given q this defines set of values of 0^';^ where sev- 
eral vacua are degenerate. This set is independent of 
number of flavours Nf. Furthermore, we are interested 
only in those 0^';^ which minimise the term linear in mass, 
which is the second condition for ^crit • It is easy to check 
that independently of q the value 6'crit ~ tt satisfies both 
conditions. Other critical values depend on topological 
charge of the domain q. In general if q = ni /n2 then there 
are ni critical values of theta. For the value q — 0.15, as 
was fit in the domain model, we find 

^'crit = {-7r/3,7r/3, tt} (mod 27r), 

where the Dashen phenomenon occurs. We conclude that 
for any Nf and rational q there is a finite number of 
critical points in the interval [0, 27r] including 9 — t:. 

IV. SUMMARY AND DISCUSSION 

The central result of this paper is that 27r-periodicity 
of amplitudes in 9 and vacuum doubling with sponta- 
neous CP violation at certain critical values of 9, in 
particular 9 = n, are obtained for any number of light 
flavours Nf and arbitrary rational topological charge q. 
This is achieved in a model whose input is a particular 
class of nonperturturbative gluon configurations. From 
this model have been derived both confinement, the cor- 
rect pattern of chiral symmetry breaking, certain static 



7 



characteristics of the vacuum (string constant, conden- 
sates, topological susceptibility) as well as properties of 
the meson spectrum (correct splitting between masses 
of iso- vector and iso-singlet states). The most impor- 
tant qualitative feature of this class of vacuum fields is 
that it can be seen as an ensemble of densely packed 
lumps of (anti-)self-dual gluon fields, characterised by a 
finite correlation length. In this paper and in previous 
works it has been shown on a semi-quantitative level that 
such vacuum fields can reproduce all of the qualitatively 
important features of the QCD vacuum associated with 
confinement and chiral symmetry realisation. 




FIG. 5: The scalar quark condensate as a function of 6 for 
Nf — 3 and q = in units of K, which is equivalent to 

corresponding plot from . The meaning of dashed and solid 
lines is the same as in Fig. 

We discuss now our results in light of the two main 
approaches, that of anomalous Ward identities and of the 
effective chiral Lagrangian at large Nc, which originally 
predicted the Dashen phenomenon in the presence of a 
theta term. For convenience, we have summarised the 
salient features of these approaches in Appendix B. 

The first point of comparison is the minimisation with 
respect to non-singlet angles. In both and 'G], this is 
referred to as "Dashen's theorem" . In the domain model 
context, precisely the same conditions Ea. H22l23|) have 
arisen by considering the free energy as a function of the 
dynamical variables /?" associated with the fiavour non- 
singlet angle dependence of fermion boundary conditions. 
In [3 and [g the corresponding angles are denoted (pi cor- 
responding to the phases of an Nj x Nf matrix, either 
the matrix of scalar condensates (tpiipj) in |^ or the chiral 
field of 0. In all cases, these results (for the degenerate 
mass case) emerge because of the structure of the sym- 
metry breaking term — toH X^tfi cos$i. The significance 
of this structure for the domain model is that it emerges 
precisely from the spectral asymmetry of the Dirac op- 
erator in the domain field background as computed via 
zeta function regularisation ■ The analogues of the 
in 1^, 1^ are our functions $i associated with the flavour 
dependence of the fermion domain boundary conditions. 

The second point of contact relates to the existence or 
otherwise of an axial U{1) Goldstone boson in the spec- 



trum. In both 1^, 1^, the input that there be no such 
boson leads to a constraint on the sum of the aforemen- 
tioned angles and the theta parameter: 



Nf 



(24) 



In the domain model approach, as mentioned at the out- 
set, the anticipated absence of such a boson is a conse- 
quence of the lack of continuous axial U{1) degeneracy of 
the ground state of the free energy for the massless case. 
But this ground state is determined by 



Nf 



arctan(tan <I>,j) = 



(25) 



from Ea. H19|) . Resolving the arctantan structure brings 
us to the same form as in Ea. H24|l . when and $i are 
identified, again consistent with our observations above. 
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FIG. 6: The pseudo-scalar quark condensate as a function of 
Q for N f — 3 and q — 1/Nf in units of H. The meaning of 
dashed and solid lines is the same as in Fig. 

Certainly the reason for these coincidences is most 
transparent in the comparison to the large Nc approach, 
reviewed in appendix B: the effective singlet and non- 
singlet meson fields there are contained in U{Nf) ma- 
trices but diagonahsed by SU{Nf) x SU{Nf) transfor- 
mations. In the domain model, where the conditions on 
quark fields at domain boundaries also involve the U{Nf) 
fiavour matrices exp[i(a-|-/3"T°/2)75] which are similarly 
diagonahsed by special unitary matrices. 

The main difference with Q is the specific value of 
the topological charge. For a more direct comparison 
with 5] we set q = l/Nf in our formula, Ea. H18|) . The 
Z— dependent term is then 27rl which is inconsequential 
under the cosine for the condensate. Thus our results 
completely agree with for this case. For the Nf — 3 
case, we plot the scalar and pseudoscalar condensates in 
Figs. Inland El which are identical to corresponding plots 
in 5]. 

The significance of our result in comparison to 0| is 
that the domain model assumes the dominance of specific 
class of nonperturbative gluon configurations in the QCD 
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functional integral and thereby provides the 2tt period- 
icity in theta dependence and the existence of critical 
values of theta parameter as output simultaneous with 
the correct resolution of the axial U{1) problem. In 
the absence of a pseudoscalar U{1) boson is input in or- 
der to obtain the theta dependence, while the responsible 
vacuum structure is unknown. 

The present work can be generalised for the case of 
non-degenerate quark masses in the presence of the theta 
term for which an analysis using anomalous Ward identi- 
ties is given in Q . Because of the above-noted similarities 
with the domain model the condensate dependence on 9 
will be identical. 

Finally, we mention that the model under considera- 
tion with rational topological charges q reflects the strong 
CP problem in the usual way, but certainly cannot resolve 
it. As in several other approaches (for instance we 
notice that the free energy in the model is minimised by 
6 = 0, but there is no reason within the model to de- 
mand the minimisation of the free energy with respect 
to 9, which is an external parameter here. However, al- 
lowing irrational values of q can drastically change the 
status of the strong CP problem in the model due to the 
appearance of infinitely many critical values of 9 in any 
finite interval. We shall analyse this intriguing possibility 
in a separate publication. 
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APPENDIX A: PSEUDO-SCALAR CONDENSATE 

We treat here Nf = 1 and the integral for one domain 
for simplicity; the generalisation is straightforward. Con- 
sider 

(V;e"^T=V)a = lim w-1Z±(to,/3) 

rn^O 

where (— )+ corresponds to an (anti-)self-dual domain. 
A chiral transformation: 

^ = e-'^-'^iP' (Al) 

for a field '0 belonging to J-a, namely satisfying 



leads to the transformed field V'' satisfying 

so that ip' belongs to Ta-p- Performed in the integral, 
this chiral transformation results in 

Z±(m,/3)=e^i^(«Z±_^(m,0) (A2) 

where the phase F is fixed by /3-independence of Z^{Q, (i) 
and the result for the fermionic determinant 

Z^(0,0) = exp(±i(7 arctan[tan(Q:)]), 

which gives 

F(/3) = ±2 arctan[tan(/3)]. (A3) 
Thus, collecting all together, we get 

Z^(rn,P) = exp {zbiq arctan[tan(a)] 
+vmexp{±i{a — . 
In the presence of the 0-term this gives 

(V!e*''^=V')a = -e=^*' ('''■'=*''"[*''"(")l~'')[cos(a-/3) 
±i sin(a — 

which finally, after summing self-dual and anti-self-dual 
configuration, leads to 

{i}e'f^^'"ll})a = -e'"™"^'? (arctan[tan(a)]-e)}^ ^^^^^ _ 

Finally using this result for the complete calculation with 
N oo domains and integration over angles aj associ- 
ated with j-th domain, we get 

(V;e*/3^=0) = -^J2coii{aM{9)-f3), 

kl 

where the sum spans all 27r-inequivalent minima aki{9) 
of the free energy density, explicitly given in Eq. H18I) . 

In particular, the pseudoscalar condensate corresponds 
to /3 = 7r/2 and thus reads 

{^ij5i^) =iV/H^sin(afc,(6')) 

kl 

It should be stressed here again, that for any 9 and f3 
J2cos{aki{9)-P) = 0. 

kl 

We see that the ^-dependence of the scalar condensate 
{f3 = 0) for different minima of the free energy den- 
sity is given by cos(afc;(6')) while the dependence of 
the pseudoscalar condensate {f3 = n/2) is described by 
sm{aki{9)), as shown in Figs. |21 and ^ 

The crucial point in this derivation is that the chiral 
transformation Eq. (|A1I) changes the space of integration 
and simultaneously generates a phase F in Eq. l|A2p . This 
phase is fixed in the form of Eq. (|A3|) by substituting 
m = into Eq. (|X2l) 

Z±(0,/?) = e^i^(^)zt^(0,0), 

and taking into account the known explicit form of 
Z^_p{0,0) and Z^{0,P), where the latter does not de- 
pend on /3 by construction. 
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APPENDIX B: OTHER DERIVATIONS OF THE 
DASHEN PHENOMENON 

1. Anomalous Ward identities 

We summarise the salient aspects of and related 
works. Denote by J^s the singlet axial vector current 
renormalised gauge-invariantly. Its divergence gives the 
anomaly. Inserting this current into Green's functions 
with a composite operator consisting of a product of local 
observables Ok{xk), taking the divergence and using the 
anomaly one obtains 

dliT{0\J^,ix)l[Ok{xkm = 

k 

2Nd^T{0\K^ix)l[Ok{xkm 

k 

I k 

The quantity in the first term of the RHS is the well- 
known Chern-Simons current arising from the anomaly. 
The second term arises from commuting the divergence 
through the T-product, which generally gives a commu- 
tator of the operators Ok with the axial charge, and then 
rewriting that commutator in terms of chiralities Xk cor- 
responding to Ok which are defined via the eigenvalue- 
like relation 

[Q5,0fc] =-XfeOfc. (B2) 

The charge does not correspond to J^5 but rather to 
the gauge-dependent, conserved axial current. Despite 
the gauge dependence of Q5, chiralities of Ea. llB2p are 
renormalisation group- and gauge-invariant 5j. For a 
right handed quark field x — I. The condition for avoid- 
ing a U{1) boson is that the LHS of Ea. (jBl|) vanishes at 
zero momentum transfer. Taken between 9 vacua, one 
can rewrite the Chern-Simons current contribution (due 
to its connection to topological charge density) on the 
RHS via a derivative with respect to 9 yielding 

= [2Nfi-^ - J2 Xi]r(0| J]Ofe(xfe)|0). (B3) 
; k 

Now one chooses the local operator product qq such that 
the sum of chiralities is xi = 2. One extracts then 
the relation 

{iNf^ - l){0\uLUn\0) = (B4) 

for example. This can be trivially solved for the conden- 
sate. The next step is to recognise the condensate here 
as an element of the matrix of condensates which break 
the non-singlet chiral symmetry 

(0|fe).(to),|0) = Cl/.,. 

This matrix can be brought by chiral rotations into a 
form diag(e"^'). The real angles (j>i parametrising the 



matrix now carry the 9 dependence of the condensate. 
Inserting this into Ea. Hij4|l gives 

1=1 

whose solution is 

^0,-0 = 0. (B5) 

i=l 

This equation should be understood here as a direct con- 
sequence of the requirement that no zero mass boson cou- 
ple to the gauge- invariant axial vector current. 

The above considerations should be repeated in the 
presence of such a perturbation matrix of masses uii for 
each flavour 

= 2 ^ rrii cos (B6) 

i=l 

According to Dashen's theorem, the true vacuum is found 
by minimising a quark mass term with respect to small 
chiral rotations about this configuration, meaning shift- 
ing under the cosine by infinitesimal LOi. Minima are 
determined by the conditions 

rrii sin 1/)^ — independent of i = A 

cos(f>^ > 0. (B7) 

To give an explicit solution assume Nf ~ 2 with degen- 
erate quark masses. Then the consequence of Dashen's 
theorem Ea. (|Ij7() gives sin0i = A/to. Solutions consis- 
tent with Eg. ljByp are 0^ = e [— 7r/2, 7r/2] modulo 27r, 
with (/) = arcsinA/m. The absence of a Goldstone bo- 
son Ea. (|B5|) gives 01-1-02 = 9. Combining these gives 
the result that the symmetry breaking term and thus the 
condensate is | cos(6'/2)|. Essentially the absolute value 
appears because the cosine may not change sign, due to 
Ea. (|B7l) . as 9 varies. Out of this emerges that the peri- 
odicity of 9 is 27r. For general iV/ with degenerate quark 
masses the corresponding result is that the condensate is 
proportional to cos[0(mod27r)/A^/]. 



2. Large Nc approach 

Now we briefly summarise how identical results are ob- 
tained in the large approach of Witten and Veneziano- 
Di Vecchia |6|. An TV/ x N f unitary field C/, parametrised 
as [/ = ?7o exp(i7r°t°) with iV| meson (non-singlet and 
singlet) fields tt", is considered: the U{Nf) generators t°- 
include the identity matrix as well as the usual SU{Nf) 
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generators. An effective Lagrangian for U with chiral and 
axial U{\) symmetry broken is 



C 



where the first term is a kinetic term, the second a mass, 
and the third term is designed to only yield a term 
quadratic in the singlet field and is consistent with large 
Nc counting rules. The mass matrix M. is then diago- 
nalised by an SU{Nf) x SU{Nf) transformation (corre- 
sponding to the action on the quark mass matrix) to the 
form 

with M = diag(/if). The are combinations of the 
squares of the meson masses but are linear in the cor- 
responding quark masses rm, oc mi, generalising 
the Gell-Mann-Oakes-Renner relationship and contain- 
ing the quark condensates - see [2^ for explicit formu- 
lae. The diagonalisation leads to a corresponding phase 
transformation on the f7— field: U e~^^/'^W such that 
the In dot term undergoes a shift 0. The aim is now to 
minimise the effective potential to find the vacuum con- 
figurations. This is aided by considering a diagonal U 
parametrised as diag(e*"^') where the (j)i are complicated 
functions of the meson fields tt". This leads to the po- 
tential 



a 



The terms with /i^ evidently arise from the separate 
phases in U and thus after expanding the (pi in pow- 
ers of the meson fields tt" will give mass terms for the 
N'j — 1 non-singlet mesons. The last term incorporates 



only information about the overall phase of U; the corre- 
sponding expansion to second order in tt"^ will reflect the 
mass of the singlet f/(l) state. That this is taken to be 
large compared to the other N'j — 1 is now input at this 
point. One now minimises with respect to the angles (pi 
("Dashen's theorem") giving 



^1 sin ( 




Evidently then one has again 



/Zj sin pi = independent of i 



identical to Ea. (|B7(l . On the other hand, taking the first 
N'j — 1 mesons to be light, the only way to minimise the 
potential for a heavy U{1) boson is to constrain the term 
with factor a/Nc to be exactly zero, namely 



= 



which is identical to the condition from anomalous Ward 
identities - no surprise since the effective Lagrangian is 
engineered to satisfy these identities. Implementing these 
conditions leads to the identical dependence on 9 dis- 
cussed in the previous section, which now appears in the 
potential at its minima 

V(0i)min ^ -F^^flf COS pi\min- 



Since the quark condensates reside in the fif an identical 
dependence on 9 as that found using anomalous Ward 
identities emerges. The dependence involves cos{9/Nf) 
now because of the way the [/—field transforms un- 
der chiral transformations in terms of the 9 parameter: 
Indet U Indet U — i9 and there is no relative factor of 
two between the sum of the phases of U and 9. 
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